Loss of superfluidity in a Bose-Einstein 
condensate on an optical lattice via a novel 
classical phase transition 

Sadhan K. Adhikari 

Instituto de Fisica Teorica, Universidade Estadual Paulista, 
01.405-900 Sao Paulo, Sao Paulo, Brazil 

Abstract 

We predict the loss of superfluidity in a Bose-Einstein condensate (BEC) trapped 
in a combined optical and axially-symmetric harmonic potentials during a resonant 
collective excitation initiated by a periodic modulation of the atomic scattering 
length a, when the modulation frequency equals twice the radial trapping frequency 
or multiples thereof. This classical dynamical transition is marked by a loss of 
superfluidity in the BEC and a subsequent destruction of the interference pattern 
upon free expansion. Suggestion for future experiment is made. 
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The detailed study of quantum phase effects on a macroscopic scale such as 
interference of matter waves [1] has been possible after the experimental load- 
ing of a cigar-shaped Bose-Einstein condensate (BEC) in a combined axially- 
symmetric harmonic plus a standing-wave optical lattice potential traps in 
both one [2,3] and three [4] dimensions. There have been several theoretical 
studies on different aspects of a BEC in a one- [5,6] as well as three-dimensional 
[7] optical-lattice potentials. The phase coherence among different optical- 
lattice sites of a trapped BEC on an optical lattice has been established in 
recent experiments [2,3,4,8,9] through the formation of distinct interference 
pattern when the traps are removed. This long-range phase coherence in the 
condensate along the entire optical lattice is a sign of communication among 
various sites which is necessary for developing superfluidity in the condensate. 

The phase- coherent BEC on the optical lattice is a superfluid [4,10] as the 
atoms in it move freely from one optical site to another by quantum tunnel- 
ing through the high optical potential barriers. It has been demonstrated in 
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a recent experiment by Greiner et al. [4] that, as the optical potential traps 
are made much too higher, the quantum tunneling of atoms from one opti- 
cal site to another is stopped resulting in a loss of superfluidity in the BEC. 
Consequently, no interference pattern is formed upon free expansion of such 
a BEC which is termed a Mott insulator state [4,10], in which an individ- 
ual atom is attached to a fixed optical site and its free mobility to a nearby 
site by tunneling is stopped as in an insulator. This phenomenon represents 
a superfluid-insulator quantum phase transition and cannot be properly ac- 
counted for in a mean-field model based on the Gross-Pitaevskii (GP) equation 
[11] where these quantum effects are mostly lost. 

Following a suggestion by Smerzi et al. [12], Cataliotti et al. [13,14] have 
demonstrated in a novel experiment the loss of superfluidity in a BEC trapped 
in a one-dimensional optical-lattice and harmonic potentials when the cen- 
ter of the harmonic potential is suddenly displaced along the optical lattice 
through a distance larger than a critical value. Then a modulational instabil- 
ity takes place in the BEC. Consequently, it cannot reorganize itself quickly 
enough and the phase coherence and superfluidity of the BEC are lost. The 
loss of superfluidity is manifested in the destruction of the interference pattern 
upon free expansion. Distinct from the quantum phase transition observed by 
Greiner et al. [4], this modulational instability responsible for the superfluid- 
insulator transition is classical dynamical in nature [12,13,15]. This process is 
also different from the Landau dissipation mechanism [12,15], occurring when 
the fluid velocity is greater than local speed of sound. The present dynam- 
ical phase transition can be well described by the mean-field GP equation 
[6,12,15]. When Landau instability occurs, the system lowers energy by emit- 
ting phonons [15]. The GP equation can not simulate energy dissipation and 
hence the Landau dissipation mechanism. 

The above modulational instability is not the unique dynamical classical pro- 
cess leading to a superfluid-insulator transition. Many other classical processes 
leading to a rapid movement or collective excitation in the condensate may 
lead to such a transition. The movement should be rapid enough so that the 
BEC cannot reorganize itself to evolve through orderly phase coherent states. 
Here we suggest that a collective resonant excitation of the BEC may also 
lead to the destruction of superfluidity due to a classical superfluid-insulator 
transition. There have been theoretical [16] and experimental [17] studies of 
collective excitation in the BEC in the absence of an optical trap initiated by a 
modulation of the trapping frequency. The study of such collective excitation 
in the presence of an optical trap has just began [18]. 

In the present study the collective excitation is initiated near a Feshbach res- 
onance [19] by a periodic modulation of the repulsive atomic scattering length 
a (> 0) via a — * a + Asin(Qt) where t is time, A an amplitude, and Q the 
frequency of modulation. Such modulation of the scattering length can be re- 
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alized experimentally near a Feshbach resonance by manipulating an external 
background magnetic field. Although, the background magnetic field and the 
scattering length are nonlinearly related in general, for small modulations A 
(< a) an approximate linear relation between the background magnetic field 
and the scattering length may hold which might make the implementation 
of the above modulation experimentally possible. When Q = 2u or multiples 
thereof, resonant collective oscillation can be generated in the BEC, where 
uj is the radial trapping frequency [20]. This resonant oscillation destroys the 
superfluidity of the BEC provided that the condensate is allowed to experi- 
ence this oscillation for a certain interval of time called hold time. We base 
the present study on the numerical solution of the time-dependent mean-field 
axially-symmetric GP equation [11] in the presence of a combined harmonic 
and optical potential traps. This transition involving collective excitation of 
the BEC and described by the mean-field GP equation is classical, rather than 
quantum, in nature. 

The time-dependent BEC wave function ^(r;r) at position r and time r is 
described by the following mean-field nonlinear GP equation [11] 
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where m is the mass and iV the number of atoms in the condensate, g = 
<iiih 2 a/m the strength of interatomic interaction. In the presence of the com- 
bined axially-symmetric and optical lattice traps V(r) = ^mw 2 (r 2 + v 2 z 2 ) + 
V^pt where uj is the angular frequency of the harmonic trap in the radial di- 
rection r, voj that in the axial direction z, with v the aspect ratio, and V^p^ 
is the optical lattice trap introduced later. The normalization condition is 
Jdr\^( r ;r)\ 2 = l. 

In the axially-symmetric configuration, the wave function can be written as 
\l/(r, r) = ip(r,z,r). Now transforming to dimensionless variables x = \/2r/l, 



y = \/2z/l, t = tuj, I = (moj), and tp(x, y; t) = x\J I 3 / \^8ip(r, z; r), Eq. (1) 
becomes [21] 
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<f(x,y;t) = 0, 



(2) 



where nonlinearity n = Na/l. In terms of the one-dimensional probability 
P(y, t) = 27r / °° dx\ip(x, y, t)\ 2 /x, the normalization is given by Jf^ dyP(y, t) = 
1. 
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We use the parameters of the experiment of Cataliotti et al. [8] with repulsive 
87 Rb atoms where the radial trap frequency was uo = 2tt x 92 Hz. The optical 
potential created with the standing-wave laser field of wavelength A = 795 nm 
is given by V ^ = V E R cos 2 (k L z), with E R = h 2 k 2 L / '(2m), k L = 2ir/\, and V 
(< 12) the strength. For the mass m = 1.441 x 10~ 25 kg of 87 Rb the harmonic 
oscillator length I = sJh/{rmJ) = 1.126 jum and and the dimensionless time 
unit u)- 1 = 1/(2tt x 92) s = 1.73 ms. In terms of the dimensionless laser wave 
length Ao = V2X/1 — 1 and the dimensionless energy Er/(TiuS) = 47r 2 /Ag, 
V opt of Eq. (2) is V ovt /(hu) = V^/Xl) cos 2 (27ry/A ). 

We solve Eq. (2) numerically using a split-step time-iteration method with 
the Crank-Nicholson discretization scheme described recently [21]. The time 
iteration is started with the known harmonic oscillator solution of Eq. (2) 
with n = 0: cp(x,y) = [^/ (87T 3 )] 1 / 4 xe~^ x +uy ^ 4 [21]. The nonlinearity n as 
well as the optical lattice strength V are slowly increased by equal amounts in 
lOOOOn steps of time iteration until the desired nonlinearity and optical lattice 
potential are attained. Then, without changing any parameters, the solution 
so obtained is iterated 50 000 times so that a stable solution is obtained 
independent of the initial input and time and space steps. 

The one-dimensional pattern of BEC on the optical lattice for a specific nonlin- 
earity and the interference pattern upon its free expansion have been recently 
studied using the numerical solution of Eq. (2) [6]. Here we study the destruc- 
tion of this interference pattern after the application of a periodic modulation 
of the scattering length resulting in a similar modulation of nonlinearity n in 
Eq. (2) via 

n -> n + Asin(Qt), (3) 



where A is an amplitude. In the present model study we employ nonlinearity 
n — 5, the axial trap parameter v = 0.5, and the optical lattice strength 
Vq — 6 throughout. First we calculate the ground-state wave function in the 
combined harmonic and optical lattice potentials. 

When the condensate is released from the combined trap, a matter-wave in- 
terference pattern is formed in a few milliseconds as described in Ref . [6] . The 
atom cloud released from one lattice site expand, and overlap and interfere 
with atom clouds from neighboring sites to form the robust interference pat- 
tern due to phase coherence. The pattern consists of a central peak and two 
symmetrically spaced peaks, each containing about 10% of total number of 
atoms, moving in opposite directions [8,9,6,22]. 

If we introduce the modulation of nonlinearity (3) after the formation of the 
BEC in the combined trap, the condensate will be out of equilibrium and 
start to oscillate. As the height of the potential barriers in the optical lattice 
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Fig. 1. One-dimensional probability P(y, t) vs. y and t for the BEC on optical lattice 
under the action of modulation (3) with n = 5, Vt = 2lo and A = 3 and upon the 
removal of the combined traps after hold times (a) 17 ms, (b) 35 ms, and (c) 52 ms. 

is much larger than the energy of the system, the atoms in the condensate will 
move by tunneling through the potential barriers. This fluctuating transfer of 
Rb atoms across the potential barriers is due to Josephson effect in a neutral 
quantum liquid [8] . The phase coherence between different optical sites of the 
condensate may be destroyed during this rapid oscillation with large amplitude 
and no matter-wave interference pattern will be formed after the removal of 
the joint traps. 

However, resonant oscillation arising from modulation (3) is not the only clas- 
sical dynamical mechanism for the destruction of superfluidity. It can also be 
destroyed via a periodic modulation of the radial trapping potential [23] or 
by giving the BEC a large displacement along the optical trap [9] . For a weak 
modulation or small displacement, the superfluidity is preserved independent 
of hold time in the modulated trap. For a stronger modulation and longer hold 



5 




Fig. 2. One-dimensional probability P(y, t) vs. y and t for the BEC on optical lattice 
under the action of modulation (3) with (a) = u and (b) £1 = 3u> with n = 5 and 
A = 3 upon the removal of the combined traps after hold time 69 ms. 

time, there is destruction of superfluidity via a classical dynamical process. 

Now we explicitly study the destruction of superfluidity in the condensate 
upon application of the modulation (3) leading to a resonant oscillation. The 
loss of superfluidity only takes place if the BEC is allowed to experience the 
resonant oscillation for a certain interval of time (hold time). The resonant 
oscillation is excited for Q = 2u or multiples thereof. Strictly speaking, this 
should happen when Q equals the natural frequency of oscillation of the system 
or multiples thereof. This natural frequency is 2uo for a noninteracting system 
with n = 0. In the presence of nonlinear mean-field interaction this frequency 
is expected to change. This change is found to be small in numerical simulation 
in the present context of small nonlinearity. However, some deviation is found 
to take place for large nonlinearity [21]. 

In this calculation we take n = 5 in Eq. (3) and for the effective nonlinearity 
after modulation to remain positive (repulsive condensate) we restrict to A < 
5. Negative values of nonlinearity corresponding to atomic attraction may 
lead to collapse and instability [11] and will not be considered here. We shall 
present results with A = 3 in this study, although any other A, which is not 
negligibly small, leads to similar result. For Q = 2uj the BEC executes rapid 
oscillation exciting collective resonant modes [17,18] resulting is a destruction 
of superfluidity. 
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For numerical simulation we allow the BEC to evolve on a lattice with x < 
20 (Ltm and 20 /im > y > —20 //m after the modulation (3) is applied and 
study the system after different hold times. The one- dimensional probability 
P(y,t) is plotted in figures 1 (a), (b) and (c) for hold times 17 ms, 35 ms 
and 52 ms, respectively. For hold time 17 ms, prominent interference pattern 
is formed upon free expansion. In figure 1 (a) three separate pieces in the 
interference pattern corresponding to three distinct trails can be identified. 
The interference pattern is slowly destroyed at increased hold times as we can 
see in figures 1 (b) and (c). As the hold time increases the maxima of the 
interference pattern mix up upon free expansion and finally for the hold time 
of 52 ms the interference pattern is completely destroyed as we find in figure 1 
(c). As the BEC is allowed to evolve for a substantial interval of time after the 
application of the periodic modulation in the scattering length, a dynamical 
instability of classical nature sets in which destroys the superfluidity [12,13]. 

The superfluidity reappears rapidly as frequency Q of the modulation (3) is 
changed to a nonresonant value. We demonstrate this for Q = uo and 3u in 
the following. In figures 2 (a) and (b) we present the evolution of probability 
P(y,t) after the application of the modulation for n = 5, A = 3 and Q — u 
and 3u. We see from figures 2 (a) and (b) that in both cases the interference 
pattern is obtained after a hold time of 69 ms. 




Fig. 3. One-dimensional probability P(y, t) vs. y and t for the BEC on optical lattice 
under the action of modulation (3) with £1 = 4lo with n = 5 and A = 3 upon the 
removal of the combined traps after hold time 69 ms. 

The resonant collective oscillation also appears for higher multiples of Q = 2uj 
thus leading to a destruction of superfluidity. This is illustrated in Fig. 3 
where we present the evolution of probability P(y,i) after the application of 
the modulation with n = 5, A = 3 and Q = 4u for a hold time of 69 ms. 
We see that no interference pattern is formed in this case upon release of the 
BEC from the joint traps. In this case the superfluidity is maintained for a 
hold time of 52 ms whereas for Q = 2uj it was destroyed at 52 ms as can be 
seen in figure 1 (c). 

In the absence of the optical trap, a cigar-shaped BEC can be excited to 
collective resonant states for modulation (3) when fl is an even multiple of uj 
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or vuj [20]. The present study shows that in the presence of an optical trap, 
only the former possibility leads to a loss of superfluidity, at least for the small 
nonlinearity (n = 5) considered here. The possibility Q = 2voj does not seem 
to lead to prominent collective resonant excitation and hence does not easily 
lead to a loss of superfluidity as one can see from figure 2 (a) with v = 0.5. 
The reason for the absence of collective resonant excitation in this case is not 
clear. It is possible that for large nonlinearity n and amplitude A in (3) there 
is breakdown of superfluidity in this case. 

In conclusion, we have studied the destruction of superfluidity in a cigar- 
shaped BEC loaded in a combined harmonic and optical lattice traps upon 
the application of a periodic modulation of the scattering length leading to 
resonant collective excitation when the frequency of modulation equals twice 
the radial trapping frequency or multiples thereof. In the absence of modula- 
tion, the formation of the interference pattern upon the removal of the com- 
bined traps clearly demonstrates the phase coherence [7,6]. At these resonance 
frequencies the phase coherence is destroyed signaling a superfluid-insulator 
classical phase transition, provided that the BEC is kept in the modulated 
trap for a certain hold time. Consequently, after release from the combined 
trap no interference pattern is formed. The superfluidity in the BEC is quickly 
restored when the frequency of modulation of the scattering length is changed 
to a nonresonant value away from Q = 2oo or multiples. It is possible to study 
this novel superfluid-insulator classical phase transition experimentally and a 
comparison of those results with mean-field models will enhance our under- 
standing of matter wave BEC. After the recent experiments of Cataliotti et 
al. [13] and Muller et al. [22] that study seems possible in a not too distant 
future. 
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